Implicit FEM-FCT algorithms and discrete Newton methods for transient convection problems
نویسنده
چکیده
A new generalization of the flux-corrected transport (FCT) methodology to implicit finite element discretizations is proposed. The underlying high-order scheme is supposed to be unconditionally stable and produce time-accurate solutions to evolutionary convection problems. Its nonoscillatory low-order counterpart is constructed by means of mass lumping followed by elimination of negative off-diagonal entries from the discrete transport operator. The raw antidiffusive fluxes, which represent the difference between the highand low-order schemes, are updated and limited within an outer fixedpoint iteration. The upper bound for the magnitude of each antidiffusive flux is evaluated using a single sweep of the multidimensional FCT limiter at the first outer iteration. This semi-implicit limiting strategy makes it possible to enforce the positivity constraint in a very robust and efficient manner. Moreover, the computation of an intermediate low-order solution can be avoided. The nonlinear algebraic systems are solved either by a standard defect correction scheme or by means of a discrete Newton approach whereby the approximate Jacobian matrix is assembled edge-by-edge. Numerical examples are presented for two-dimensional benchmark problems discretized by the standard Galerkin FEM combined with the Crank-Nicolson time-stepping.
منابع مشابه
A semi-implicit FEM-FCT algorithm for efficient treatment of time-dependent problems
A new generalization of the flux-corrected transport (FCT) methodology to implicit finite element discretizations is proposed. The underlying high-order scheme is supposed to be unconditionally stable and produce time-accurate solutions to evolutionary convection problems. Its nonoscillatory low-order counterpart is constructed by means of mass lumping followed by elimination of negative off-di...
متن کاملAlgebraic Flux Correction I Scalar Conservation Laws
This chapter is concerned with the design of high-resolution finite element schemes satisfying the discrete maximum principle. The presented algebraic flux correction paradigm is a generalization of the flux-corrected transport (FCT) methodology. Given the standard Galerkin discretization of a scalar transport equation, we decompose the antidiffusive part of the discrete operator into numerical...
متن کاملExplicit and implicit FEM-FCT algorithms with flux linearization
A new approach to the design of flux-corrected transport (FCT) algorithms for linear/bilinear finite element approximations of convection-dominated transport problems is pursued. The raw antidiffusive fluxes are linearized about an intermediate solution computed by a positivity-preserving low-order scheme. By virtue of this linearization, the costly evaluation of correction factors needs to be ...
متن کاملMassively Parallel Solution Techniques for Higher-order Finite-element Discretizations in CFD
The purpose of this paper is to present techniques to solve higher-order finite element discretizations on massively parallel architectures. Implicit schemes are considered as a means of achieving mesh independent convergence rates for both time dependent problems and steady state solutions obtained through pseudo-transient continuation. Domain decomposition preconditioners are presented for th...
متن کاملOn the design of flux limiters for finite element discretizations with a consistent mass matrix
The algebraic flux correction (AFC) paradigm is extended to finite element discretizations with a consistent mass matrix. A nonoscillatory low-order scheme is constructed by resorting to mass lumping and conservative elimination of negative off-diagonal coefficients from the discrete transport operator. In order to recover the high accuracy of the original Galerkin scheme, a limited amount of c...
متن کامل